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1.  Background 

Composite  materials  are  being  increasingly  used  throughout  the 
Army  as  a  means  of  decreasing  the  weight  and  increasing  the 
strength  of  Army  systems.  With  the  increased  use  of  composite 
materials,  the  question  of  electromagnetic  shielding  must  be 
answered.  The  Army  Research  Laboratory  (ARL)  is  currently 
involved  in  a  composite  shielding  program  whose  objective  is  to 
examine  the  composites  used  in  Army  systems  and  characterize 
them  electromagnetically.  The  approach  used  has  been  to  develop  a 
composite  characterization  laboratory  that  will  measure  the  electri¬ 
cal  properties  of  composites.  The  Helmholtz  coil  is  one  tool  that  will 
be  used  in  the  characterization  laboratory. 


2.  Introduction 

Helmholtz  coils  are  used  to  develop  uniform  magnetic  fields  for  a 
variety  of  research  applications.  Helmholtz  coils  can  be  easily  con¬ 
structed  and  the  fields  easily  calculated.  This  makes  them  especially 
useful  in  calibrating  sensors  and  numerous  other  low-frequency 
magnetic-field  testing  applications.  This  report  discusses  the  calcula¬ 
tion  and  measurement  of  Helmholtz  coil  fields  to  be  used  for  testing 
of  composite  materials  to  low-frequency  magnetic  fields.  Future 
composite  testing  of  this  nature  will  be  based  on  American  National 
Standard  ANSI/ASTM  A698-74,  "Standard  Test  Method  for  Mag¬ 
netic  Shield  Efficiency  in  Attenuating  Alternating  Magnetic  Fields." 


3.  Test  Configuration 

The  Helmholtz  coils  wfere  constructed  in  accordance  with  the  guide¬ 
lines  in  the  ASTM  698-74  standard.  They  were  fabricated  from  3/»-in. 
plywood,  V:-in.  Extren®  threaded  rod,  and  14  AWG  copper  (coated) 
wire.  Each  coil  had  a  diameter  of  approximately  47.5  in.  (1.219  m), 
with  48  turns  each  of  the  copper  wire.  The  coils  are  shown  in  figure 
1.  Sixteen  V:-in.  Extren  threaded  rods  were  used  to  form  the  circular 
shape  of  the  coil.  The  magnetic  fields  associated  with  the  Helmholtz 
coils  were  measured  with  two  single-turn  shielded  loops.  Both  were 
constructed  from  0.141-in.  copper  jacketed  cable,  one  with  a  diam¬ 
eter  of  2  in.  and  the  other  of  3  in.  Additional  information  on  the  cali¬ 
bration  of  the  loops  is  provided  in  section  4. 

A  Hewlett  Packard  (HP)  3577A  network  analyzer  was  used  to  meas¬ 
ure  the  magnetic  fields  of  the  Helmholtz  coil.  The  output  of  the  ana¬ 
lyzer  was  connected  to  an  audio  amplifier,  whicn  was  connected  to 
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the  Helmholtz  coil.  The  reference  channel  monitored  the  current  in 
the  loop.  Channel  A  recorded  the  voltage  associated  with  the  2-  or  3- 
in.  loops  placed  inside  the  coil.  The  complete  test  configuration  is 
shown  in  figure  2. 


Figure  1.  Helmholtz 
coil  construction. 


Figure  2.  Test 
configuration. 


Helmholtz  coils 
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4.  Measurement  Loop  Calibration 

Two  electrostatically  shielded  single-turn  loops  (one  2-in.,  one  3-in.) 
were  used  to  measure  the  fields  associated  with  the  Helmholtz  coils. 
Both  were  constructed  as  shown  in  figure  3.  The  open-circuit  voltage 
at  the  loop  output  terminals  can  easily  be  calculated  by 

Vx  =  2 nfNAB  ,  (1) 

where 

Voc  =  open-circuit  voltage  at  loop's  output  terminals, 

/  =  frequency  (in  hertz), 

N  =  number  of  turns  in  loop  (one  turn), 

A  =  area  of  loop  (in  meters),  and 
B  =  magnetic  flux  density  (in  teslas). 

The  loops  were  calibrated  in  a  transverse  electromagnetic  (TEM)  cell 
known  as  a  Crawford  cell.  Using  an  HP  35 77 A  network  analyzer,  we 
measured  the  loop  voltage  as  a  function  of  output  voltage  applied  to 
the  TEM  cell  versus  frequency.  Knowing  the  physical  dimensions  of 
the  TEM  cell  and  the  measured  area  of  the  loops,  we  could  deter¬ 
mine  the  ratio  of  the  measured  versus  reference  voltages;  to  deter¬ 
mine  the  calculated  response,  we  applied  this  ratio  to 

Rc  =  201og(Vv/Prt,)  ,  (2) 


Figure  3.  Measurement 
loop  construction. 
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where 


Rc  =  calibrated  response, 

V  =  open-circuit  voltage  at  loop's  output  terminals,  and 
Vrgf  =  reference  voltage  applied  to  the  TEM  cell. 

The  measured  and  calculated  results  are  shown  in  table  1  for  100 
kHz.  Figure  4  compares  the  measured  and  calculated  responses  for  a 
frequency  range  of  100  kHz  to  10  MHz  for  both  the  2-  and  3-in. 
loops.  One  can  see  that  the  calculated  and  measured  responses  fol¬ 
low  each  other  to  well  within  0.5  dB.  Based  on  equation  (1),  it  is 
assumed  that  both  of  the  loops  would  continue  to  fall  at  20  dB  per 
decade  to  dc.  Because  of  limitations  in  the  measurement  system,  the 
loops  could  not  be  calibrated  below  100  kHz. 


5.  Helmholtz  Coil  Predicted  Results 

The  magnetic  fields  generated  by  a  circular  wire  loop  can  be  derived 
according  to  the  following  procedure.  The  electromagnetic  (EM) 
fields  generated  by  a  circular  wire  loop  carrying  current  /  will  satisfy 
the  Maxwell  equations.  For  a  wire  loop  (or  a  coil)  excited  with  low- 
frequency  current,  almost  all  the  energy  is  stored  in  the  magnetic 
field.  This  energy  is  determined  from  the  concepts  and  rules  of  cir¬ 
cuit  theory.  Circuit  theory  can  be  regarded  as  describing  a  special 
class  of  solutions  of  the  Maxwell  equations  that  result  from  using  the 


Table  1.  Loop 
calibration  results. 


Loop 

Measured 

Calculated 

TEM  cell  measured 

dimension 

area 

response 

response 

Difference 

(in  ) 

(m) 

(dB) 

(dB) 

(dB) 

2 

2.545 

-95.1511 

-94.9368 

0.2143 

3 

5.572 

-88.3447 

-87.8906 

0.4541 

Figure  4.  Measured  (solid)  and  calculated  (dotted)  responses  of  (a)  2-  and  (b)  3-in.  loops. 
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first  term  of  a  power  series  solution  for  the  fields.  The  fields 
computed  with  this  approximation  are  called  quasi-static.  EM 
phenomena  such  as  radiation  are  neglected  in  the  quasi-static 
approximation.1 

A  general  rule  of  thumb  for  circuit  elements  excited  by  sinusoidal 
voltage  or  current  is  that  the  quasi-static  approximation  will  be 
extremely  good  when  the  physical  size  of  the  element  is  sufficiently 
small  with  respect  to  the  wavelength  of  the  source  excitation.  This 
restriction  will  be  satisfied  if  the  phase  of  the  voltage  (current)  is 
approximately  the  same  over  the  spatial  extent  of  the  element.  For  a 
circular  loop,  this  condition  requires  that  the  current  be  nearly  con¬ 
stant  over  the  circumference  of  the  loop.  The  current  can  change  in 
time,  but  at  each  instant  of  time,  the  current  in  every  part  of  the  loop 
should  have  the  same  value.  Therefore,  the  quasi-static  approxima¬ 
tion  will  be  very  good  as  long  as  2na  «  A  (or  in  terms  of  frequency 
instead  of  wavelength,  v  «  c/lOa),  where  a  is  the  loop  radius  and  c 
is  the  velocity  of  light  in  free  space. 

One  can  find  the  EM  fields  generated  by  a  loop  in  this  order  of 
approximation  by  determining  the  current  in  the  loop  employing 
circuit  theory  and  using  that  uniform  current  as  the  source  of  the 
vector  potential  from  which  the  fields  can  easily  be  derived.  This  is 
the  approach  followed  below.2 

Figure  5  shows  the  coordinates  and  variables  for  a  loop  of  radius  a 
located  in  the  xy  plane.  The  loop  is  assumed  to  have  a  total  current  I. 
The  vector  potential  for  a  general  current  distribution  can  then  be 
written  as 


where  A  is  the  vector  potential  and  /  is  the  current  in  the  loop.  From 
the  direction  of  the  current,  A  has  only  a  0  component,  and  by  sym¬ 
metry  A  is  independent  of  the  variable  0.  In  setting  up  the  integra¬ 
tion,  we  arbitrarily  set  0  =  0  to  simplify  the  results.  With  the  symbols 
given  in  figure  5,  the  following  can  be  derived: 


dl  =  (-a  sin  0',  a  cos  0\  0 )d0'  , 
r  =  (r  sin  6,  0,  r  cos  6)  , 
r  =  [a  cos  0',  a  sin  0\  0)  , 

|  r  -  r  |  =  >J rl  +  a1  -  2 ra  sin  6  cos  0'  . 


M.  Fano,  L.  J.  Chu,  and  R.  B.  Adler,  Electromagnetic  Fields,  Energy,  and  Forces,  John  Wiley  &  Sons,  Inc.  (1960). 
2William  R.  Smythe,  Static  and  Dynamic  Electricity,  Hemisphere  Publishing  Corp.  (1989). 
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Figure  S.  Coordinates 
and  variables 
describing  a  circular 
loop  used  to  express 
EM  fields. 


Substituting  these  variables  in  the  integral,  we  note  that  the  integral 
involving  sin  0'  vanishes  (it  is  odd),  as  it  should  since  only  an  A $ 
component  must  survive,  and  the  integral  involving  cos  0'  (even) 
allows  the  range  of  integration  to  be  halved.  This  results  in 

AJir.8>=Ugl  f  . .  . 

2n  J0  s/r2  +  a1  -  2 ra  sin  6  cos  0' 

It  is  more  convenient  to  express  this  integral  in  terms  of  cylindrical 
coordinates  using  r2  =  p2  +  z2  and  sin  6  =  p/<J p 2  +  z 1  ,  which  gives 

A  (07)  =  H»!A  f  cos  0'd0' 

*  2ft  Jo  Jp2+z2  +  aZ_  2 ap  COS  0' 

This  integral  does  not  have  a  closed  form;  however,  there  is  a  trans¬ 
formation  that  results  in  tabulated  functions.  Change  the  variable  to 
make  the  upper  limit  of  the  integral  n/2,  by  using  0’  -  ft  +  20.  The 
integral  becomes 


P _ L2_!i£l 

0V/VJ  -  n  r - - - 7 

Jo  v(a+P)  +z* 


-  1  )d0 


-  Aap  sin:  0 


If  we  define  k 2  = 


Aap 


(< a+pf  +  z 2 


,  the  integral  is  transformed  into 


10 


-w*) ■  ^  4%  [( 1  -  2  k’)K(k>  -  £H 


with 


W)  = 


d0 


o  l -k:  sinT 


and  E(k) 


k *  sin'  0  d0  . 


K  is  the  complete  elliptic  integral  of  the  first  kind,  and  £  is  the  com¬ 
plete  elliptic  integral  of  the  second  kind.  Both  are  tabulated 
functions. 

The  magnetic  fields  are  computed  from  the  vector  potential  by 
B  =  V  x  A,  giving  two  components: 


BP(P,z)  =  -~  and  B:(pJ  =  j,Jt(pA0)  . 

To  calculate  these  derivatives,  we  will  need  the  derivatives  of  K,  E, 
and  k,  since  they  are  functions  of  p  and  z.  The  math  is  straightfoward 
and  the  results  are  summarized  below: 


3 K  E 

K 

„„h  M  = 

3 it  k(l-k2) 

k 

dk 

k  k 

dk  _  k  k 3 

k 3 

dp  2  p  4  p 

4a 

dz 

4a  p 

The  final  expression  for  the  magnetic  field  components  becomes 


P\J  {P  +  a)2  +  z2 
1 


lK  \J{p  +  af  +z2 


y-P'  +  \E(k)-K(k) 
(' a-py  +  z 2 

a2-p2-z 2 


(a-p)2  +  z2 


E(k)  -  K(k) 


For  a  single  loop  in  the  xy  plane,  we  can  find  the  magnetic  field  com¬ 
ponents  at  any  location  p  and  z  by  first  determining  k  from  the  above 
definition,  and  then  substituting  into  the  above  formulas. 

We  can  check  the  accuracy  of  these  expressions  for  the  simple  case  of 
the  magnetic  field  components  along  the  axis  of  a  loop.  On  the  axis, 
p  =  0,  which  implies  k  =  0,  and  hence  fC(0)  =  E(0)  =  n/2.  The  expres¬ 
sion  for  Bp  becomes  indeterminate,  but  application  of  l'Hopital's 

rule  gives  Bp  =  0.  The  expression  for  Bz  becomes  0.5/V?://  y  {a2  +  z2f , 
as  expected. 

Calculating  the  magnetic  fields  from  a  Helmholtz  coil  requires  fields 
computed  for  two  loops  having  the  same  axis,  located  in  the  planes 
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at  z  =  ~d  and  z  =  +d.  The  separation  distance  2d  is  related  to  the  loop 
radius  by  d  =  a/ 2.  The  expression  for  the  total  magnetic  field  written 
in  terms  of  the  above  derivation  is  just  the  sum  of  the  individual 
fields  from  each  loop.  Using  the  notation  and  results  derived  earlier 
for  an  individual  loop,  one  can  write 

B';'al(p,z)  =  Bp(p -  +  d)  +  BP(P~  -  d)  , 

BT*{pz)  =  B.(p,z  +  d)  +  Bz(p  -  -  d)  . 

Once  the  solutions  were  obtained,  the  variables  were  substituted  for 
the  actual  Helmholtz  coil.  The  diameter  of  the  coil  was  measured  at 
1.219  m,  the  average  current  measured  (see  sect.  6)  was  0.019052  A, 
and  the  number  of  turns  was  48. 

6.  Measured  Results 

Numerous  magnetic  field  measurements  were  taken  on  the  test  con¬ 
figuration  described  in  section  3.  The  objective  of  the  magnetic  field 
measurements  was  to  verify  the  predicted  response  and  determine 
the  size  of  the  test  volume.  This  was  done  by  essentially  mapping 
the  inside  of  the  Helmholtz  coil.  A  cylindrical  coordinate  system  was 
used  for  consistency  with  the  predicted  results.  Both  the  p  and  z 
coordinates  were  varied  in  2-in.  increments.  Because  of  symmetry, 
the  <j>  component  was  not  recorded,  and  only  one  quadrant  was 
measured. 

Figure  6  is  a  measurement  of  the  Helmholtz  coil's  fields  as  a  function 
of  frequency.  Figure  6  clearly  shows  that  the  coil  has  a  resonance  at 
approximately  58  kHz,  which  would  mean  that  the  usable  frequency 
range  of  the  Helmholtz  coil  is  between  dc  and  40  kHz.  Frequencies 
above  40  kHz  should  not  be  used. 

The  measured  results  compared  very  well  to  the  calculated  results, 
as  is  shown  in  figures  7  to  9.  Figure  7  shows  a  three-dimensional 
comparison  between  the  measured  and  calculated  z-component  of 
the  magnetic  field  as  a  function  of  p  and  z.  One  can  see  that  the  two 
compare  very  well  in  wave  shape  and  magnitude.  Figure  8  shows  a 
two-dimensional  comparison  of  the  measured  versus  calculated 
magnetic  field  as  a  function  of  z  on  the  centerline  (p  =  0). 

The  peak  magnitude  of  the  measured  field  is  1.345  x  10-6  T  (teslas), 
while  the  calculated  is  1.348  x  10-6  T.  Both  are  essentially  the  same 
magnitude.  The  measured  data  have  a  noticeable  peak  at  z  =  0, 
where  the  calculated  data  are  flat  out  to  z  =  ±0.1.  This  difference 
(0.37  percent)  cannot  be  fully  explained,  but  one  possibility  is  that 
the  measurement  loop  was  not  perpendicular  to  the  incident  mag- 
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Figure  6.  Measured 
resonance  of 
Helmholtz  coil. 


netic  field.  Equation  (1)  assumes  that  the  magnetic  field  is  normal  to 
the  plane  of  the  measurement  loop.  When  the  measurement  loop  is 
not  normal  to  the  incident  field,  equation  (1)  has  to  be  multiplied  by 
the  sine  of  the  angle  between  the  plane  of  the  loop  and  the  incident 
field.  A  difference  of  0.37  percent  is  equivalent  to  the  2-in.  loop  being 
off  less  than  5°  or  %2  in.  above  and  below  the  normal  plane.  While 
great  care  was  taken  to  map  the  fields  inside  the  coils,  the  measure¬ 
ment  loops  were  only  sighted  in  for  being  perpendicular,  and  errors 
of  this  magnitude  were  possible.  Figure  9  is  another  representation 
of  the  same  data,  shown  as  contour  plots.  The  contour  plots  show 
that  the  field  is  uniform  at  1.3  x  10~*  T  for  most  of  the  inner  region. 
Again,  one  can  see  that  the  measured  data  compare  very  well  with 
the  calculated. 
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Figure  9.  Contour 
plots  of  B(:l  versus  : 
and  fr.  (a)  measured 
and  <b)  calculated. 
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7.  Conclusions 


The  objective  of  this  report  is  to  provide  a  technical  guideline  for  the 
use  of  the  Helmholtz  coil  that  was  designed  and  built  for  the  ARL 
electromagnetic  composite  testing  laboratory.  The  experimental  ver¬ 
sus  theoretical  results  were  compared  for  the  primary  field  compo¬ 
nent  and  proved  to  show  excellent  correlation.  By  measuring  the 
current  being  applied  to  the  coils  and  using  the  equations  in  section 
5,  the  user  can  predict  the  magnetic  field  at  any  point  inside  the 
coils.  To  assist  the  user  in  determining  the  magnetic  fields,  we  have 
included  appendix  A,  which  is  a  Matlab®  program  listing  that  can 
be  used  to  generate  the  graphs  as  seen  in  section  6. 

The  size  of  the  test  volume  is  also  an  important  com  'ation  when 
working  with  a  Helmholtz  coil.  The  ANSI/  ASTM  A.  4  standard 

indicates  that  the  test  volume  is  related  to  the  coil  .meter.  The 
standard  states,  “The  Helmholtz  coil  diameter  shall  be  at  least  three 
times  the  length  of  the  test  specimen  or  four  times  its  diameter,  or 
both."3  Figure  9  confirms  this  test  volume,  in  that  the  field  is  uni¬ 
form  to  just  over  ±0.21  m  on  the  centerline  of  the  p-axis  (0.34  percent 
of  the  diameter)  for  all  values  of  z.  However,  figure  9  also  sh<  vs  that 
the  test  volume  depends  on  the  geometry  of  the  test  object.  For 
objects  that  are  less  than  0.28  m  high  (z-axis)  and  on  the  centerline, 
the  test  area  extends  to  ±0.4  m  (p-axis).  Objects  that  are  more  than 
0.28  m  high  have  a  test  area  of  ±0.21  m  (p-axis).  As  a  rule  of  thumb, 
the  coil  diameter  should  be  at  least  three  times  the  length  of  the  test 
specimen.  The  test  area  can  be  expanded  out  to  about  ±0.35  m 
(p-axis)  on  centerline  for  all  values  of  z,  if  a  2-dB  change  in  field  uni¬ 
formity  can  be  tolerated  at  the  comers  of  the  test  volume.  The  reader 
is  referred  elsewhere4  for  a  complete  discussion  of  how  one  can 
extend  test  volume  beyond  the  above  limits  for  the  testing  of  large 
objects. 

The  loading  of  the  coils  by  the  test  object  has  not  been  addressed  in 
this  paper.  The  user  must  be  aware  that  a  large  object  of  magnetic 
material  placed  inside  the  coil  will  load  the  coil  and  concentrate  the 
fields  on  the  test  object.  The  resultant  field  will  be  distorted  when 
compared  to  the  field  when  the  coil  was  empty.  It  is  recommended 
that  the  user  measure  the  current  in  the  empty  coil  and  verify  that  it 
does  not  change  when  the  test  object  is  inserted.  If  the  current  does 
change,  the  amplifier  must  be  adjusted  to  the  original  level. 


3ANSI/ASTM  A698-74,  American  National  Standard,  Standard  Test  Method  for  Magnetic  Shield  Efficiency  in 
Attenuating  Alternating  Magnetic  Fields  (1974). 

4£.  L.  Bronaugh,  Helmholtz  Coils  for  EMI  Immunity  Testing:  Stretching  the  Uniform  Field  Area,  Seventh  Interna¬ 
tional  Conference  on  Electromagnetic  Compatibility  (Conf.  Publ.  No.  326). 
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Appendix  A.  Matlab  Program  Listing  for  Field 
Calculations 

The  following  is  a  Matlab  program  used  to  calculate  the  fields  at  any 
point  inside  a  Helmholtz  coil  of  Radius  a  and  current  /. 

%  Like  loops. a  but  symmetric  plan*  through  z  axis 
%  computes  the  B  field  from  2  thin  loops  of  rsdius  a, 

%  locatsd  at  s  -  +d  and  -d,  (d-a/2,  Bslaholtz  Coils),  and 
%  currant  X. 

LM  —  48  %  Number  of  turns  in  the  loops 

auo  -  4*pi*le-7> 

I-. 019052; 
a-. 6096; 
d-a/2; 

%  Define  the  plane  over  wnich  fields  are  computed 
%  M  aust  be  odd  to  include  the  point  (0,0) 

M— 26;  %  No.  of  points  along  the  rho  axis 

N-51;  %  No.  of  points  along  the  z  axis 

Pl-linspace(0,a,M) ;  %  rho  always  positive 

P  -  [  PI (M:-l :2 )  PI  ]; 

Z-linspace ( -d , d , N ) ; 

[p,z]  -  aeshgrid ( PI , Z ) ;  «  rho-0,  p(:,l);  z-  -d,  z ( 1 , : ) 

«  Determine  modulus  and  Elliptic  Integrals 
kl-(4*a*p) ./(  (a+p) . *2  ♦  (z-d). “2  ); 

[K1,B1]  -  ellipke(kl); 

k2-(4*a*p) ./(  (a+p). *2  ♦  (z+d).“2  ); 

[K2.E2]  -  ellipke(k2) ; 

%  Compute  B-rho  fields 

bpl  -  ( (z-d) ./(p.*sqrt(  (a+p). *2  +  (z-d). “2  )))  .*  ... 

(  (  (a“2  +  p. *2  ♦  (z-d) .“2) ./( (a-p) .“2  +  (z-d). '2  )  ).*E1  -  K1  ); 
bp2  -  ( (z+d) ./(p.*sqrt(  (a+p). *2  ♦  (z+d).*2  )))  .+  ... 

(  (  (a*2  +  p.*2  ♦  (z+d) .*2) ./( (a-p) .*2  +  (z+d). “2  )  ).*B2  -  K2  ) ; 
bp  -  (  (LN»muo*I)/(2*pi)  )* (bpl  +  bp2); 

bp(l:N)  -  zeros(N,l);  %  Remove  NaN's  from  rho-0  axis 

bp  -  (  bp(t,M:-ls2)  bp  J ;  %  Make  it  symmetric 

%  Compute  B-z  fields 

bzl  -  (  1  ./sqrt(  (a+p). *2  +  (z-d). *2  )  )  .*  ... 

(  (  (aA2  -  p.*2  -  (z-d)  .*2) ./( (a-p)  .*2  +  (z-d). “2  )  ).*E1  +  K1  ); 
bz2  -  (  1  ./sqrt(  (a+p). *2  +  (z+d).A2  )  )  .*  ... 

(  (  (a“2  -  p.*2  -  (z+d). *2). /((a-p). *2  +  (z+d). *2  )  ).*E2  +  K2  ); 

bz  -  (  (LN*muo*I)/(2*pi)  )*(bzl  +  bz2); 

bz  -  [  bz(:,M:-l:2)  bz  ];  %  Make  it  symmetric 

%  Compute  the  total  B  field 
bt  -  sqrt (bp. *2  +  bz.*2); 

%  test  along  the  z  axis  where  we  know  the  solution 
bzz-  1  . /(sqrt(a“2  +  (Z-d) . A2) ) .“3  +  1  ./(sqrt(aA2  +  (Z+d) . A2) ) .*3; 
bzz  -  0.5*LN*muo*I*aA2*bzz; 
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